Parity-time symmetric optical coupler with birefringent arms 
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In this work, we propose a parity-time iVT-) symmetric optical coupler whose arms are birefrin- 
gent waveguides as a reahstic physical model which leads to a so-called quadrimer i.e., a four complex 
field setting. We seek stationary solutions of the resulting linear and nonlinear model, identifying 
its linear point of VT symmetry breaking and examining the corresponding nonlinear solutions that 
persist up to this point, as well as, so-called, ghost states that bifurcate from them. We obtain 
the relevant symmetry breaking bifurcations (between symmetric -circularly polarized- and asym- 
metric -elliptically polarized- states) and numerically follow the associated dynamics which give 
rise to growth/decay even within the "PT-symmetric phase. Our symmetric stationary nonlinear 
solutions are found to terminate in saddle-center bifurcations which are analogous to the linear VT- 
phase transition. We found that the VT symmetry significantly changes the stability and dynamical 
properties of the modes with different polarizations. 

PACS numbers: 42.65. Jx, 42.65.Tg, 42.65.Wi 
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I. INTRODUCTION 

An optical coupler with one arm having losses and an- 
other one having gain, balanced against each other, re- 
cently became a test-bed for many phenomena originat- 
ing from the interplay of the parity-time (VT) symme- 
try and nonlinear! ty. Unidirectional dynamics un- 
versality of the dynamics 0, symmetry breaking proper- 
ties d, [1] , switching of the beams Q and of solitons @ , 
formation of symmetric and asymmetric bright solitary 
waves [3js|, breathers 0], and their stability dark 
solitons [11 1, as well as the emergence of ghost states [l2l- 
[T3 | and large-scale temporal T'T-symmetric lattices [l5| 
are some among the many topics that have been touched 
upon in very intense recent theoretical and experimental 
work. 

As a direct extension of the previous activity, a large 
chunk of which has focused on the prototypical setting 
of the PT-symmetric dimer, there emerges a problem of 
effect of VT symmetry and nonlinearity on the polariza- 
tion of the electric field. In that regard, the previously 
proposed settings, to the best of our knowledge, were 
chiefly focused on effectively scalar models. On the other 
hand, the "vector" type of problems is natural for exper- 
imental settings where the exploited fibers obey birefrin- 
gence, since the two orthogonal polarizations are to be 
taken into account [l6| . 

In this work we consider a T^T-symmetric coupler 
whose arms are birefringent waveguides. Assuming that 
the first waveguide is active and the second one is ab- 
sorbing, we address the problem of a PT-preserving (in 
the linear limit) configuration. While being an interest- 
ing model from a physical point of view, this setting also 
offers a different (in comparison to what was studied be- 
fore) mathematical situation where the nonlinear modes 
bifurcate from doubly degenerated eigenvalues of the lin- 
ear problem. This requires the generalization of earlier 



developed approaches (e.g. like the one reported in [TtI ) 
for the bifurcation of the nonlinear modes from the lin- 
ear spectrum. In addition, it presents a rich playground 
for dynamical systems analysis, due to the emergence of 
a variety of saddle-center bifurcations (nonlinear analogs 
of the linear T'T-phase transition), as well as symmetry- 
breaking (pitchfork) ones. It is these nonlinear states, 
their emergence, stability, dynamics and the asymptotics 
of the system that we will focus on hereafter. 

The organization of the paper is as follows. In sec- 
tion ini we present the model in its evolution as well as 
in its stationary form. In section IIIIl we focus on its 
linear properties. Then, in section IIVI we examine the 
nonlinear modes (and bifurcations). Finally, in sectionlVl 
we briefly touch upon the dynamical implications of our 
findings and in section IVll we present some conclusions, 
as well some potential directions for future work. 



II. THE MODEL 

We specify the problem by imposing that the principal 
optical axes of the two Kerr- type waveguides are 7r/4- 
rotated with respect to each other, as it is schematically 
represented in Fig. [T] In each arm, labeled by j, there 
arc two orthogonal field components of the electric fields 
which we write down in the form [l6| (j = 1,2): 



Ej(r,2;,t) = [uj{z)Aj{v -Tj 
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+ u,+,{Mo+2{v - r,)e-'ft+^^e,+2] J -e^^' + c.c. (1) 
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Here Uj are the field envelopes depending on the prop- 
agation distance z, i.e. we consider the stationary - in 
time - problem, assuming that the carrier wavelength 
Ao is in the region of the normal group velocity disper- 
sion, thus ruling out a possibility of modulational insta- 
bility; r = [x, y) is a transverse radius vector, and ri.2 
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are the positions of the centers of the cores of the cou- 
pler. The real parameters /3j are the propagation con- 
stants of each of the field components, and are the 
polarization vectors, which are mutually orthogonal in 
each arm of the coupler, i.e. ei • 63 = 62 ■ 64 = 0. The 
real functions Aj{r — rj) and Aj+2(r — r^) describe the 
transverse distributions of the fields in each waveguide 
and the normalization coefficient where x is the 

Kerr coefficient, is introduced for convenience. For j ^ I 
and j = 2 the functions Aj{r) are centered in different 
points Vj. Also, for the sake of simplicity, we consider 
Aj{r) = Aj+2(r) = A{r) (for j = 1,2), such that the 
integral 



J A{r~r,)A{r-r,+i)d\ 
/A2(r)d2r 



(the integration is performed over the transverse plane) 
describes the linear coupling between the respective 
modes. Since in the configuration shown in Fig. [1] 
ei • 62 = ei • 64 = 63 • 64 = —63 • 62 = I/V2 we use the 
single linear coupling coefficient k (see also [l3)- Then 
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FIG. 1: (Color online) (a) Schematic presentation of a VT- 
symmetric coupler based on birefringent fibers, (b) Equiv- 
alent graph (plaquette) representation illustrating the VT- 
symmetry. Here — and -I- stand for active and lossy waveg- 
uides, respectively. 

following the analysis described in details in [l^ we end 
up with the system of equations: 



.dui 
dz 



-fc(ii2 + U4) + - ( |wi| +-|ti3| ]ui 

-lulule^"^^^ (2a) 

dU2 A l2 2 I |2\ 

= -fc(ui - U3) - 171*2 - MM2I +g|W4| I U2 
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w^<e-*^^" (2d) 
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Here 7 > describes gain in the first waveguide and dis- 
sipation in the second waveguide, Aj = (/3^ — P'j+2) 



with /3j = Wo being the carrier wave frequency, 

is a properly normalized mismatch between the propaga- 
tion constants of the orthogonal polarizations Uj+2 and 
Uj. The asterisk stands for complex conjugation. 

We will be interested in the stationary solutions, in 
particular in their linear stability properties and ensuing 
nonlinear dynamics which can be found in the two proto- 
typical limiting cases of (i) zero mismatches — /3j+2| = 
and (ii) large mismatches — li'jj^2 \ ^ fcAo/c when the 
respective nonlinear terms can be neglected. It is conve- 
nient to introduce a parameter a which vanishes {a = 0) 
in the case (ii) and is unity (a ~ 1) in the case (i). Using 
the standing wave ansatz Uj{z) — Wje'^^^, where Wj are 
2;— independent, into ([2]), we obtain the system of alge- 
braic equations: 



bwi = k{w2 + W4) - i7U'i + ( + -1^3^ ) wi 



a 



wlwl, (3a) 



bw2 = k{wi - W3) + i"fW2 + ( |lf2p + 3l'"^4p ) W2 



a 



w^wi,, (3b) 



bws = k{wi - W2) - i-/W3 + (^gkiP + k3pj W3 

+ ^wlwl (3c) 

bW4 = k{wi + W3) + ijWA + ^-|W2P + |w4p^ ^4 

+ ^wlwl. (3d) 

Below the spectral parameter b will be also referred to as 
the propagation constant. 



III. PROPERTIES OF THE LINEAR PROBLEM 

First we address the underlying linear problem [which 
corresponds to the situation when all cubic terms in ^ 
are negligible]. It can be rewritten in the matrix form 
few ~ Hvf where 

-i"f k k 
k 17 — fc 
— fc —17 k 
k k i"f / 

(hereafter we use tilde in order to distinguish eigenvalues 
and eigenvectors of the linear problem). 

The operator H is VT symmetric, which means that 
[H, VT] = HVT-VTH = 0, where T' is a spatial rever- 
sal linear operator 





/Wl\ 




w — 


W2 


and H = 1 
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and T performs element-wise complex conjugation: 
Tw = w* . The spectrum of operator H consists of two 
double eigenvalues 



(5) 



which are real for 7 < 7cr'' where 7^"* = \f^k will be re- 
ferred to as a primary critical point: the spontaneous VT 
symmetry breaking occurs at 7^^ above which the eigen- 
values are all imaginary. In order to visualize the VT 



symmetry of the linear system, following 17[ one can rep- 
resent it with a graph shown in the right panel of Fig. [1] 
reminiscent of four linearly coupled waveguides [l7| (no- 
tice however the sign difference in the coupling constants) 
or plaquettes [T9| . 

Since the details of our analysis are the same for both 
eigenvalues 6_ and 6+, we drop the subscripts + and 
— wherever this does not lead to confusion. Despite 
having double eigenvalues in its spectrum, H is diago- 
nalizable below the T^T-symmetry breaking point. This 
means that double eigenvalues are scmisimple, i.e. for 
an eigenvalue b one can find two linearly independent 
eigenvectors, i.e. H-w^^'> = b\v^^\ where j = 1,2. More- 
over, each eigenvalue b possesses an invariant subspace 
spanned by w*^^) and -w^^y 

Let us also notice the following peculiarity of the case 
at hand. In a situation where a T^T-symmetric operator 
H has no multiple eigenvalues, the condition of unbroken 
VT symmetry (i.e. reality of all the eigenvalues) requires 
that for each eigenvalue b the corresponding eigenvector 
w can be chosen as an eigenstate of the PT operator, 
i.e. VTyv = w. However, in the situation at hand arbi- 
trarily chosen linearly independent eigenvectors w'-'^-' and 
w'2) may not be VT cigenstatcs. However, unbroken VT 
symmetry requires that a certain linear combination of 
w'^^ and jg an eigenstate for the VT operator. More 
specifically, it is easy to establish that all the eigenvec- 
tors that belong to the invariant subspace of b and, at the 
same time, arc the cigenstatcs for the VT operator, can 
be parametrized by a complex parameter a as follows: 



k _ " 



V 



(6) 



Being interested in linearly independent vectors w, it is 
sufficient to consider only the vectors with 

a = e'^. (7) 

Then Eq. ([6|) yields a monoparametric set of eigenvec- 
tors •w{6) with a real parameter 9. In particular, setting 

6'^^^ = and 6^"^^ ~ arctan(^4p^) one can choose two or- 
thogonal (and therefore linearly independent) eigenvec- 
tors: 



(hereafter we use the standard scalar product (g,h) = 
Yl'^j=i93^*j)- eigenvector w corresponding to the 

eigenvalue b can be represented as a linear combination of 
VT eigenstates w'-^'^): w = Aiw'^-* + A2w''2) (of course, 
this does not mean that any eigenvector w is also a VT 
eigenstate). 

Let us also introduce a Hermitian adjoint operator H'^ . 
Since the matrix H is symmetric, one has H"^ = H* . As 
long as VT symmetry of H is unbroken, the spectrum 
of the adjoint operator H'^ also consists of two double 
eigenvalues b± which are semi-simple. Any eigenvector 
corresponding to an eigenvalue b of the adjoint operator 
can be represented as a linear combination of (w^^^)* 
and (w(2))*. 



IV. NONLINEAR MODES 

A. Bifurcations from the linear limit 

Now we develop a perturbation theory for the eigen- 
states of the linear problem giving rise to monoparamet- 
ric families of nonlinear modes. We will look for nonlin- 
ear modes w that are eigenstates of the VT operator, i.e. 
VT'w = w. To this end we introduce the expansions 



w = e\v{9) + e-'^Wg 



and b = b 



(9) 



Here e is a small real parameter, W3 and B2 are the 
coefficients of the expansions and 6 is to be determined 
from the symmetry of the solution (see below) . We notice 
that for the expansion to be meaningful the coefficient B2 
must be real. 

Expansions © describe nonlinear modes that bifur- 
cate from the linear limit corresponds to e = and is 
given by the eigenvector 'w{0) being a linear combination 
of such that VT'w{9) = w(0). Respectively, in the 

linear limit the parameter b is given by the eigenvalue b. 

Passing from e = OtoO<e^l one has to compute 
the coefficients W3 and B2. While the physical sense 
of the coefficient B2 is clear — it is a deviation of the 
propagation constant due to small nonlinearity, it turns 
out that B2 also has a clear geometrical interpretation. 
Indeed, let us consider the total energy flow through the 
coupler, which is defined by 



U 



E 



(10) 



Expansion ([9]) implies that in the vicinity of the linear 
limit U =^ e2(w(0),w(e')) + 0{e^). Therefore, the co- 
efficient B2 governs a slope of the energy curve in the 
vicinity of the bifurcation point, i.e. 



W(1'2) 



(w(l),w(2)) 







(8) 
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(11) 
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For the sake of definiteness, now we concentrate on 
the case a = 0. Then the nonhnear problem ([3]) is con- 
veniently written in the matrix form 



bw = Hvf + _F(w)w, 



(12) 



where F(\v) is a diagonal matrix-function describing 
the nonlinearity: F{w) = diagdwij 



W2 



12 2| 



Izwal', ||?«2r + l'"^4rj- Substituting 



into Eq. ([H]), noticing that F(w) = e^F{w) + 0{e^), and 
collecting the terms order of e"^, we obtain 



{H - fe)W3 = - [F{w{e)) - B2]W(6l). 



(13) 



Equation (|13p implies two possibilities. The first one 
corresponds to the case when at some 9 the eigenvector 
w(6') of the operator H is simultaneously an eigenvector 
for the matrix F{'w{9)). Then the coefficient B2 can 
be chosen as an eigenvalue of F{'w{9)) corresponding to 
the eigenvector w{6) (provided that this eigenvalue is 
real). In this situation, the right hand side of Eq. (fT3|) 
is zero and it is sufficient to set W3 — 0. Since the 
matrix F{'w{9)) is diagonal, its eigenvalues are equal to 
its diagonal elements and the corresponding eigenvectors 
are given as columns of the 4x4 identity matrix. 

Let us first assume that all the eigenvalues of F{'w{9)) 
are simple. In this case \v{9) can not be an eigenvec- 
tor for F{w{9)). The latter fact becomes evident if one 
notices that w(0) has no zero entries for any 9 [see the 
definition Therefore, w{9) can be an eigenvalue for 
F(w(0)) only if F(w(0)) has a multiple eigenvalue. Then 
one could search for a linear combination of eigenvectors 
corresponding to the multiple eigenvalue, which would be 
equal to w{9). However, using the same argument, i.e. 
the fact that all entries of the vector w{9) are nonzero, 
one can see that even if F{'w{9)) has a double or an triple 
eigenvalue, the matrix F{w{9)) still can not have \v{9) 
among its eigenvalues. Therefore w(6') can be an eigen- 
vector of F{'w{9)) only if all its eigenvalues are equal. 
Imposing this constraint on the matrix _F(w(0)), one ob- 
tains \wi \ = Iwsl and |u'2| = 1104!. Noticing that the form 
ofw(0), i.e. Eqs. © and ([7]), implies that |ti;i| ~ \w4\ — 1 
we conclude that \v{9) is an eigenvalue of F{w{9)) only 
if the moduli of all the entries of w(0) are equal to unity. 
Then matrix F{w{9)) is equal to the 4x4 identity ma- 
trix multiplied by 5/3. Requiring the moduli of all the 
entries of ■w{9) to be equal, one arrives at the equation 
for 9 whose root is given as 



. TT 1 b 

9 = arctan — 

8 2 



(14) 



Since now the moduli of all the entries of \v{9) are equal 
to unity and therefore {w{9),w{9)) = 4, i?2 = 5/3 and 
Eq. ([TT|) readily yields that in the vicinity of b the slope 
5C^/5^lb=&.c/=o = 12/5 = 2.4. Notice that the found 
value does not depend on k or 7, and thus these modes 
correspond to their counterpart in pure conservative cou- 
pler with birefringcnt arms with 7 = [Eq. (jl4[) is valid 
in this case since arctan(±cx)) = ±7r/2]. 



Let us now consider the second possibility to fulfill 
Eq. p^ . If for some 9 the corresponding w{9) is not an 
eigenvector for F{\v{9)), then one must satisfy Eq. 
choosing nonzero W3. Then the coefficient B2 is to be 
determined from the solvability condition which requires 
the right hand side of Eq. ([T^ to be orthogonal to all the 
eigenvectors of the invariant subspace of b in the spec- 
trum of the adjoint operator W . As we have established 
in Sec. IIIIl any eigenvector of from the invariant sub- 
space of b can be represented as a linear combination of 
(w^^))* and (w^^))*. Requiring the right hand side of 
Eq. p3p to be orthogonal to an arbitrary linear combi- 
nation of (w^^^)* and (w'^^)*, we arrive at the following 
relations: 

(F(w(e'))w(e'),(w(i))*) _ (F(w(e'))w(e'),(w(2))*>) 



(w(6l),(w(i))* 



Ke),(w(2))* 



= ^2. 

(15) 

In Eq. (jl5p the first equality sign is an equation which 
is to be solved with respect to 9. Once a root 9 of the 
latter equation is found, then B2 is given from the second 
equality sign. Notice that despite the fact that the vector 
■w{9) is complex, the coefficient B2 will be real [T7| . 

Substituting the expression for w{9) into Eq. (|15p . one 
obtains a rather cumbersome equation, which, however 
can be attacked with a computer algebra program. Af- 
ter some transformations, Eq. ([T5|) yields the following 
condition: 



2ie\ 



I7- V-^^ +7^1 



(16) 



The latter equation has a real root 9 only if < 7 < fc, 
i.e. when the expression under the radical is not positive. 
This result suggests that there exists a critical value of 
the gain-loss parameter which we term as the secondary 

f2) 

critical point 7cr — k, such that for sufficiently small 

(2) 

7, namely, < 7 < 7cr , there exists another family bi- 
furcating from the eigenvalue b of the linear spectrum. 

(2) 

However, this family disappears for 7 > 7cr , in spite of 
the fact that the VT symmetry of the underlying linear 
problem remains unbroken, i.e. 7cr^ < 7cr^. It is impor- 
tant to point that the newly found family of solutions 
does not correspond to equal amplitude among the dif- 
ferent nodes, and hence pertains to an elliptically (rather 
than circularly) polarized family of modes. 

Computing the corresponding value of B2 for the newly 
found elliptically polarized family, one finds that the co- 
efficient B2 does depend on k and 7 (in contrast to the 
above considered circularly polarized family, character- 
ized by i?2 =5/3 for any combination of k and 7 which 
does not violate VT symmetry). In particular, when 7 
approaches k, the coefficient B2 of the elliptically polar- 
ized family tends to 5/3, which suggests that at 7 = fc 
the circularly polarized and elliptically polarized families 
globally merge. 
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B. Algebraic analysis and numerical results 

1. Exact solutions 

Having explored nonlinear modes close to the linear 
limit, where the amplitudes of the modes are small and 
therefore they can be analyzed by means of perturbation 
theory, let us now consider nonlinear modes of arbitrary 
amplitudes (turning again to the general case a ^ 0). 
Relying on results of the previous subsection we firstly 
search for nonlinear modes which have equal intensities in 
all four waveguides. Making the substitution W2 = —iwl 
and imposing the condition Wi = w^, W3 = 11:2, which is 
necessary for a nonlinear mode to be an eigenstate of VT 
thus leading to the circularly polarized light in each of 
the coupler arms, system ^ yields the single (complex) 
algebraic equation 



a = 0, 0<7<A; 



a=l, 0<7</c 



bwi = k{l — i)w* + 



5 — a , 



Iwil Wi — IJWi. 



(17) 



Representing wi — pe"^, we obtain a bi-quadratic equa- 
tion for p yielding two families of the modes bifurcating 
from the eigenvalues b±, given by ([S|), of the linear spec- 
trum: 

3{b-h) .2.^+ _ ^±(l-0-7(l + ») 



b — a 



2k 



1 P±e''^ 








f>± 




4>± 1 



Respectively, the nonlinear modes have the following 
form: 



(19) 



Using Eqs. (|18p one can easily obtain continuous fam- 
ilies of nonlinear modes that can be identified as a func- 
tion of the propagation constant 6, for given fc and 7. 

Let us now recall that Eq. predicts that for a = 
and 7 < 7cr , there exist families of elliptically polarized 
modes having different absolute values of the polarization 
vectors. Such families were indeed found in our numerics. 
However, all such modes turned out to be unstable (see 
Figs. [5H3] and discussion below). 

In the case of zero propagation constant mismatch, i.e. 
when a = 1 , one also can find families which have differ- 
ent amplitudes of the polarization vectors. The explicit 
expressions for the families bifurcating from h± read 



wx^wX = p±e"^, W2^w\^ (-1 ± \/2)p±e 

b-b± 1 .7 

;= , (I) = T — arcsm — . 

4T2V2' 2 V2fc 



■i(p 



pI 



(20) 



Remarkably, these modes, which also describe propaga- 
tion of elliptically polarized light, are stable in a certain 
range of the parameters. 
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FIG. 2: (Color online) Prototypical examples of families of 
nonlinear modes in the plane (6, C/) for fc = 1 and gain-loss 
parameters 7: 7 = 0.5 (the upper panels), 7 = 1.1 (the lower 
panels). Left and right columns correspond to a = and 
a = 1. Stable and unstable modes are shown by solid and 
dashed lines, respectively. The families with the circular and 
elliptical polarization (if any) are marked with labels "(c/' 
and "(e/'i respectively (in the color online version families 
with the circular and elliptical polarization are also shown 
by blue and red lines, respectively). Notice that the families 
corresponding to the modes having exact solutions (|18|l and 
(|20p are represented by the straight lines. 



2. Families of nonlinear modes 

The results of our analysis of the families of nonlin- 
ear modes are summarized in Fig. [21 The upper panels 
show that for < 7 < fc each eigenvalue of the linearized 
problem gives rise to two distinct (circularly and ellipti- 
cally polarized) families of nonlinear modes. In the case 
of a = the slopes of the dependencies U (b) are close 
for the families of both types, and the elliptically polar- 
ized families are always unstable while circularly polar- 
ized families have both stable and unstable solutions. For 
a = 1 one can find stable solutions both for the families 
with circular and for those with elliptical polarization. 

For fc < 7 < v^fc, the case a = does not allow 
for elliptically polarized families [see Fig. [2]C], a feature 
which is in accordance with the perturbation approach 
developed above. In this case one can only find circularly 
polarized modes, which are unstable. On the other hand, 
for a = 1, stable and unstable modes of both types can 
be found [see Fig.[2]D]. 

Summarizing at this point, we have identified 4 sets 
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of solutions, two circularly polarized with equal ampli- 
tude at the nodes, and two elliptically polarized with 
unequal such amplitudes. These all degenerate into the 
two distinct eigenvalues b±, given by ([S]), of the linear 
problem. The circularly polarized solutions are more ro- 
bust, while the elliptically polarized ones are always un- 
stable for a ~ and stable only for small enough am- 
plitudes for a — 1. Among the circularly polarized ones, 
for < 7 < ^/2k the more fundamental state (stem- 
ming from the negative eigenvalue at the linear limit) is 
always the stable ground state of the system in contin- 
uations over the parameter 6, while the excited state is 
only stable for small enough amplitudes. 

3. Continuation over 7 

An alternative and perhaps even more telling way to 
illustrate the above features stems from fixing some value 
of 6, starting from the Hamiltonian limit of 7 = and 
subsequently identifying branches of the nonlinear modes 
by means of changing 7, as shown in Fig. [31 It is impor- 
tant to note that this alternative viewpoint affords us the 
ability to visualize bifurcations that cannot be discerned 
over variations in b. 

The relevant results for parametric continuations over 
7 are given in Figs. typical examples of the corre- 
sponding linearization spectra for different values of 7 can 
be found in Fig. [5l Here, it can be seen that a lower am- 
plitude and a higher amplitude symmetric (i.e., equal am- 
plitude) branch exist, for fixed b, from the Hamiltonian 
limit of 7 = and all the way up to the linear T'T-phase 
transition point 7^'' = V^k. At that point, the two sym- 
metric branches collide and disappear in a saddle-center 
bifurcation which can be thought of as a nonlinear ana- 
log of the linear T'T-phase transition . An additional 

very interesting feature arises precisely at the point jcr 
[see (fTS]) and the related discussion] , where it can be seen 
that both branches of equal amplitude between the sites 
become dynamically unstable for a = 0. In fact, it is seen 
that for the larger amplitude branch (associated with the 
blue circles), one pair of unstable eigenvalues arises, while 
for the smaller amplitude (red diamond) branch, two such 
pairs accompany the symmetry breaking bifurcation oc- 
curring at this critical point. A closer inspection reveals 
that the symmetric branch (blue circles) is destabilized 
through a subcritical pitchfork bifurcation with its "cor- 
responding" asymmetric state (i.e., the one degenerate 
with it in the linear limit). In the case of the lower 
amplitude (excited) state for the same 6, the situation 
appears to be more complex. In particular, there exists 
once again a subcritical pitchfork with the correspond- 
ing asymmetric branch, yet this would justify one pair 
of unstable eigenvalues and we observe two. This is be- 
cause at the same point, there also exists a supercritical 
pitchfork, which gives rise to the so-called ghost states, 
denoted by magenta plus symbols. These states are anal- 
ogous to the ones to analyzed in , but remarkably 



are not stationary states of the original problem, yet they 
are pertinent to its dynamical (instability) evolution and 
for this reason they will be examined in further detail 
separately in the dynamics section below. 

In the case of a = 1, only one pair of unstable eigen- 
values emerges for the lower amplitude branch at the 

(2) 

secondary critical point of •fcr (while the larger ampli- 
tude branch remains stable throughout the continuation 
in 7). Hence, in this case, once again a saddle-center bi- 
furcation will mark the nonlinear T'T-phase transition, 
yet the number of unstable eigendirections of each sym- 
metric branch (fundamental and excited) is decreased by 
one (0 and 1 real pairs instead of 1 and 2, respectively, for 
a = 0). In this case, in fact, both asymmetric branches 
persist up to the linear T^T-phasc transition (rather than 
terminate in a subcritical pitchfork as above), and col- 
lide and disappear with each other. Interestingly all 3 
branches (the lower amplitude, excited symmetric one 

and the two asymmetric ones) become unstable at the 
.12) 

secondary critical point jcr = k, which again points 
to the existence of corresponding ghost states. For the 
lower amplitude symmetric branch, the bifurcating ghost 
states are again identified by the magenta plus symbols 
in Fig. m 




FIG. 3: (Color on-line) The four panels denote the solu- 
tion amplitude (top left), phase differences between adjacent 
nodes (top right), real and imaginary parts (second row) of 
eigenvalues for q = 0, 6 = 2, and fc = 1. The blue circles 
branch corresponds to the "+" sign in Eq. H18|) . while the 
red diamonds branch corresponds to the "-" sign (the sym- 
metric/circularly polarized branches). The green stars and 
black squares crosses are those solutions with distinct abso- 
lute values of the polarization vectors (the asymmetric or el- 
liptically polarized branches). In the top left panels, they 
collide and disappear in two subcritical pitchfork bifurcations 
with the blue circles and red diamond branches, respectively. 
The magenta pluses branch in the panels represents the ghost 

state solutions, which bifurcate from the red diamonds at 

(2) 

7 ~ "Icr = 1 and terminate at 7 = 2.44. 



V. DYNAMICS OF THE POLARIZATION 

To examine the potentially symmetry breaking (and 
more generally instability driven) nature of the dynami- 
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FIG. 4: The four panels show the same diagnostics as in the 
previous figure but now for a — 1, 6 = 2, and k = 1. 
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FIG. 5: StabiUty plots. The top two panels are for a = 0, 
and the bottom one is when a = 1. In the case a = 0, 

(2) 

at 7 = 7cr = 1, one pair of blue circles and two pairs of red 
diamonds collide at so that one pair of real eigenvalues arises 
in blue circles branch whereas two pairs of real eigenvalues 
arise in the red diamonds branch. The asymmetric branches 
only exist (and are unstable) for the smaller value of 7 = 0.5, 
for a = 0. For the same parameters (& = 2, fc = 1), in the 
case of a = 1, the excited symmetric and both asymmetric 
branches are unstable for 7 — 1.2. 



cal evolution past the critical points identified above, we 
have also performed direct numerical simulations which 
are illustrated in Fig. [51 see also Fig. [T] Here, it can 
be seen that although the relevant parameters are be- 
low the critical point for the linear T'T-phase transi- 
tion 7^' = V^k, nevertheless, symmetry breaking phe- 
nomena are observed due to the dynamical instability 
of the relevant states (the ones denoted by dashed lines 
in Fig. [5]). This dynamics may, in principle, be associ- 
ated with the so-called ghost states of complex propaga- 
tion constant that have recently been proposed as rele- 
vant for the dynamical evolution in . To substantiate 
this claim, we note that it is observed in the left panel 
of Fig. [H] that the relative phase of the two gain sites 
that lock into an equal growing amplitude, is 7r/2, as is 



those of the decaying amplitude lossy sites. In light of 
this, we seek ghost states with precisely this phase dif- 
ference and are able to explicitly identify them via the 
ansatz ^3 = iwi, W4 = {1112, setting Wj = Cje**J for 
j = 1, 2. For these branches, the propagation constant is 
complex. This highlights the potential growth or decay 
of such states. Importantly also, note that these states 
arc "ghosts" because they may be solving the stationary 
problem of Eqs. (|3]), but the U(l) invariance of the orig- 
inal model does not permit them to be a solution of the 
dynamical Eqs. 

The algebraic conditions that this family of solutions 
satisfies are 



sm0fc = ; 2 , 2 
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(24) 

(25) 
(26) 



Notice that the imaginary part of the propagation con- 
stant Bsin{(j)b) is proportional to the difference — c^. 
Hence, prior to the symmetry breaking, the relevant so- 
lutions bear a real propagation constant. Past the bi- 
furcation point one (unstable) branch has c| > , while 
the stable branch has c\ > c\. The relevant ghost state 
branches and their bifurcation from the equal amplitude 
ones are explored in Fig. [3]|4l Given that these are only 
ghost solutions of the original dynamical problem, the 
interpretation of their linearization spectrum (shown for 
completeness in Fig. is still an open problem. 

These ghost states appear, in fact, to be the ones di- 
rectly observed in the unstable evolution dynamics. To 
illustrate this, we observed the particular behavior of the 
unstable modes and how it depends on the form of the 
initial perturbation. A typical example in which the gain 
sites lead to growth and the lossy sites to decay is shown 
in the left panel of Fig. [SI It is remarkable that the 
dynamics appears to closely follow the ghost states iden- 
tified above. This is clearly showcased in Fig. [71 through 
the comparison of the growth pattern in the relevant sites 
(and the decay pattern in the lossy sites) with the ex- 
act, shifted in the propagation distance to fit the onset 
of growth, ghost state solution for the same parameters. 
The matching is excellent, lending strong support to the 
dynamical relevance of the ghost states. On the other 
hand, in the right panel of Fig. [5] a different scenario of 
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evolution is illustrated. Instead of the gain nodes grow- 
ing and the lossy ones decaying, a breathing oscillation 
settles between the two pairs. These two scenarios, illus- 
trated in Fig.[Sl are the prototypical instability evolution 
ones that we have obtained in this system. 




FIG. 6: (Color online) Dynamics of an unstable circularly 
polarized mode at & = 3, 7 = 0.5 and for a = for two dif- 
ferent small initial perturbations. The result of the evolution 
of the left panel involves growth at the gain sites (ui,3, red 
curves in the color online version) and decay at the lossy sites 
(u2,4, blue curves in the color online version). Notice that 
intensities among the two gain sites and among the two lossy 
sites are approximately equal ~ Wsf and |'U2p ~ l'"4|^) 

and are not distinguishable in the scale of the plots. In the 
right panel only the initial stage of the found persistent pe- 
riodic dynamics is shown; the simulations were performed up 
to z = 2000. 




FIG. 7: The dynamical semi-log plot of the ghost state with 
a = 0, b — 2, and k = 1 for 7 = 1.02. The red sohd line and 
blue dashed line correspond to the red diamond and magenta 
plus branches in Fig. [S] respectively. The blue dashed line is 
plotted with a shift in time, i.e. delay by z = 19.61. The 
time axis in the plot is the actual time of the red solid line. 
Their dynamical behaviors are essentially identical (and can 
not be distinguished on the scale of the plot), showcasing how 
the unstable dynamics of the red diamond branch follows the 
bifurcating ghost state. 



VI. CONCLUSION 

In conclusion, in the present work, we have proposed 
a novel, physically realistic variant of a VT symmetric 
dimer where the effect of birefringence has been taken 
into consideration. The existence of polarization of the 
electric field within the coupler yields two complex dy- 



namical equations for each of the fibers, providing a phys- 
ical realization of a plaquette model with both linear and 
nonlinear coupling between the elements. The station- 
ary states of the model were identified and both lin- 
ear and nonlinear T'T-phase transitions were obtained. 
The degenerate nature of the linear limit complicated 
the problem in comparison to other ones studied earlier 
in this context. Furthermore, the emergence of symme- 
try breaking phenomena and associated (subcritical or 
supercritical) pitchfork bifurcations, as well as their dy- 
namical implications in leading to indefinite growth and 
decay (of the corresponding waveguide amplitudes) were 
elucidated. A connection was given to ghost states that 
appeared to dominate the symmetry-breaking dynamics. 

From the physical point of view, we emphasize that the 
use of VT symmetry significantly changes the stability of 
the modes with different polarizations. For example, in 
the case of large propagation constant mismatches, we 
saw that the symmetric -circularly polarized- states be- 
come destabilized in the presence of gain/loss and may 
even cease to exist past a certain critical point of the 
relevant parameter. Instead of them, the dynamics may 
lead to a breathing exchange of "mass" between the gain 
and the lossy waveguides, or most commonly an indef- 
inite growth of the former at the expense of the latter. 
This is in line with what is expected from the dynamics of 
the dominant emerging ghost states of the system. The 
instabilities and associated dynamics should be observ- 
able in suitable generalizations of existing experiments 
such as Q . Such properties are relevant also to possibili- 
ties of solitonic -waveguide array- generalizations of the 
coupler system as well as towards a possible use of VT 
symmetric coupler for measurement techniques based on 
the use of several nonlinear modes. 

From the mathematical point of view, we believe that 
these studies may pave the way for considering multi- 
component, as well as multi-dimensional (generalizing 
the plaquettes considered here or those of [l0|) lattice 
models of PT-symmetric form. In generalizing to multi- 
plaquette configurations, it would be especially interest- 
ing to examine which of the symmetry-breaking and non- 
linear T^T-phase transition phenomena examined herein 
are preserved and what new phenomena may arise as ad- 
ditional degrees of freedom are added. Such studies will 
be deferred to future publications. 
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